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SO.  ABSTRACT 


*ttr  "T  w*€»  1 


We  have  made  significant  developments  in  the  area  of  quantum  chemical  methodology 
The  majority  of  these  developments  were  aimed  at  providing  our  group  with  state- 
of-the-art  tools  for  use  on  the  difficult  chemical  problems  in  which  we  are 
interested^-ijaking  use  of  our  new  theoretical  tools,  we  hsye  studied  the 
electron! qJsRPueT^s'1 xs tab i  li ty ,  bonding  character,  shapes)*^of  the  ground  and 
excited  electronic  states  of  several  classes  of  molecular  ions  and  novel  reactive 
intermediates. 
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During  the  three  years  this  U. S.  Army  Research  Office  Contract,  we 
have  achieved  even  more  than  was  included  in  our  original  research  pro¬ 
posal.  The  attached  list  of  published  articles  which  have  resulted  from 
Army  support  attest  to  this  fact. 

We  have  made  significant  developments  in  the  area  of  quantum  chemi¬ 
cal  methodology.  The  majority  of  these  developments  were  aimed  at  providing 
our  group  with  state-of-the-art  tools  for  use  on  the  difficult  chemical 
problems  in  which  we  are  interested.  Drs.  Robert  Donnelly,  Ron  Shepard, 
and  Ajit  Banerjee  have  been  instrumental  in  these  achievements. 

Making  use  of  our  new  theoretical  tools,  we  have  studied  the  electronic 
studies  (stability,  bonding  character,  shapes  -  of  the  ground  and  excited 
electronic  states  of  several  classes  of  molecular  ions  and  novel  reactive 
intermediates.  Included  in  our  investigations  were 


1)  anions  of  polar  molecules  such  as  LiH  ,  BeO  ,  LiF  ,  and  NaH~.  These 
Ions  have  stable  2£  ground  states  and  metastable  2x  excited  states. 

2)  Small  cluster  anions  such  as  ,  Be2~»  ®e3  »  ®e4  »  ^83  *  ^S4_» 

and  Be  Mg”.  The  alkaline  earth  clusters  Involve  the  Interaction  of  closed- 

2 

shell  atoms  (e.g.,  Mg  is  3s  )  which  gives  rise  to  very  weak  interaction 

energies  for  the  neutral  clusters  (e.g..  Mg,,).  However,  the  addition  of 

an  "extra"  electron  to  the  lowest  empty  bonding  orbital  (e.g.,  3p  for 

© 

Mg2>  gives  rise  to  Increased  interaction  energy  and,  hence,  increased 
stability  for  the  cluster  anion. 

3)  Unsaturated  carbenes  which  possess  low-energy  singlet  and  triplet  states 
whose  chemical  reactivity  patterns  vary  substantially.  ^CC:,  HjCCC:, 


,  and  H2CCCC:  have  all  been  studied.  Our  goal  was  to 


understand  the  effect  of  delocalization  and  aromaticity  on  the  relative 
energetic  stability  of  the  singlet  and  triplet  states. 
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Line  shapes  of  charge-transfer  spectra 

Jeanne  McHale*'  and  Jack  Simons6* 

Department  of  Chemistry,  University  of  Utah,  Salt  Lake  City,  Utah  84 1 12 
(Received  9  January  1079) 

The  theoretical  description  of  condensed-media  line  shapes  of  electronic  transitions,  as  developed  in  carte 
papers,  is  applied  to  the  charge-transfer  (CT)  spectra  of  the  complexes  mesitylenenetracyanoethylene 
(MesrTCNE)  and  anthracenc:chk>ranil  (ANCA).  Results  are  compared  with  those  previously  published 
for  benzenedodine  (Bz:I,),  and  the  factors  which  contribute  to  line  broadening  (in homogeneous 
broadening,  the  Franck-Condon  profile,  and  donor-acceptor  vibrational  motion)  are  -ftt-iinni  The 
results  of  an  attempt  to  determine  the  contributions  to  the  line  shape  made  by  intermolecular  orbital 
couplings  are  presented.  In  addition,  the  Bz:l:  CT  spectrum  is  reinvestigated,  and  the  question  of  which 
ionization  potential  of  benzene  is  involved  is  addressed. 


I.  INTRODUCTION 

In  earlier  publications1-1  a  theoretical  formalism  for 
describing  electronic  absorption  spectra  in  condensed 
media  was  developed.  By  taking  a  suitable  average4 
(over  the  solvent  molecules’  coordinates)  of  the  electric 
dipole  correlation  function,  the  intensity  of  absorption 
Hu)  can  be  decomposed  into  zeroth-,  first-,  and  second- 
order  contributions,  each  of  which  contains  the  effects 
of  localized  intramolecular  dipole  transitions  as  well  as 
those  of  intermolecular  charge  transfer.  The  first-  and 
second-order  contributions,  which  modify  the  zeroth- 
order  dipole  strength,  contain  the  so-called  "resonance” 
and  "fluctuation”  Integrals,  respectively,  wtych  depend 
upon  the  coupling  between  orbitals  on  neighboring  mole¬ 
cules  and  fluctuations  thereof.  The  theory  has  been  suc¬ 
cessfully  applied,  in  a  parameterized  form,  to  the  spec¬ 
tra  of  excess  electrons  in  ethanol  glass,1  anthracene 
anions  in  glassy  anthracene,1  styrene  and  cycloocta- 
tetraene  anions,1  NOJ  in  several  strongly  interacting 
solvents,1  and  to  the  benzene  iodine  charge-transfer 
spectrum.1  In  the  charge -transfer  (CT)  case,  the  zeroth- 
order  contribution  to  the  line  shape  was  shown  to  be 
given  by 
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for  a  donor  (D)  at  the  origin  and  an  acceptor  (A)  at  R. 
Ft.u(VA,  VM  VD,  VD)  is  a  product  of  the  Franck-Condon 
factors  for  the  ionization  of  the  donor  molecule  to  give 
the  donor  cation  and  the  addition  of  an  electron  to  an 
acceptor  molecule  to  give  the  acceptor  anion.  r9l),i  „A, 
is  the  transition  dipole  between  the  donor  molecular  or¬ 
bital  6  and  the  acceptor  molecular  Orbital  o,  £®°  is  the 
initial  vibrational  energy  of  the  unionized  donor  and  ac¬ 
ceptor,  A£«.a(R)  and  A£,  are  the  electronic  and  vibra¬ 
tional  excitation  energies,  respectively,  and  G( 

<r0J  is  a  Gaussian  of  width  oOA  centered  at  u\  The  in¬ 
homogeneous  broadening  constant,  e0A,  can  be  ex¬ 
pressed  in  terms  of  fluctuations  In  the  energy  levels  of 
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the  donor  and  acceptor  molecular  orbitals  as  eras  shown 
in  Ref.  1. 

The  first-  and  second-order  parts  of  Ru)  have  been 
found  in  this  laboratory  to  make  negligible  or  small  con¬ 
tributions  to  the  line  shape  in  the  cases  studied  previ¬ 
ously.  in  order  to  further  clarify  the  importance  of 
these  terms,  a  study  of  the  concentration  dependence  of 
the  Mes :  TONE  CT  spectrum  was  undertaken.  Since  the 
donor  for  this  complex  is  totally  miscible  in  common 
inert  solvents,  and  since  the  existence  of  a  ternary 
complex  involving  two  donor  molecules  has  previously 
been  reported,1  it  was  hoped  that  by  increasing  the 
mesiiylene  concentration,  die  effects  of  coupling  between 
orbitals  on  different  donor  molecules  could  be  observed. 
In  another  attempt  to  search  for  such  donor-donor  or¬ 
bital  interactions,  the  spectrum  of  AN:CA  was  also  re¬ 
corded,  but  in  a  more  limited  concentration  range. 

It  is  important  to  note  that  Eq.  (1)  gives  the  absorp¬ 
tion  intensify  for  a  fixed  donor-acceptor  separation  R. 
However,  an  important  reason  for  the  experimental  fact 
that  typical  CT  spectra  are  very  broad  is  that,  due  to 
the  Coulomblc  interaction  term  in  the  excited  state  ener¬ 
gy,  the  absorption  frequency  is  strongly  dependent  on  die 
donor -acceptor  Intermolecular  distance.  It  has  been 
found  in  this  laboratory  that,  while  Eq.  (1)  alone  is  suf¬ 
ficient  to  explain  the  width  of  spectra  which  arisq  from 
localized  transitions  in,  for  example,  styrene  and  cyclo- 
octatetraene  anions,  It  is  unable  to  fully  account  for  the 
width  of  CT  spectra.  In  order  to  model  the  homogeneous 
broadening  brought  about  by  intermolecular  vibrational 
motion,  a  thermal  average  of  the  above  intensify  expres¬ 
sion  must  be  performed,  using  a  Boltzmann  weighting 
involving  the  ground-state  donor-acceptor  potential  de¬ 
scribed  here  in  terms  of  the  harmonic  force  constant, 
k,  for  donor-acceptor  stretching: 
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R  in  this  equation  is  the  charge  separation  in  the  excited 
state,  which  is  assumed  to  be  equal  to  the  dooor-aecep- 
tor  center-of-mass  distance  for  the  f-i*  complexee 
Mes :  TCNE  and  AN :  CA.  R^  is  the  value  of  R  corre¬ 
sponding  to  the  equilibrium  geometry  of  die  donor-ac¬ 
ceptor  ground  state,  and  0  Is  (ft*  TTl.  The  factor  R1 
results  from  Mulltken’s  approximation  that,  neglecting 
overlap,  the  electronic  transition  dipole  is  proportional 
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(Reprinted  from  the  Journal  of  the  American  Chemical  Society,  101, 6174  (1979).] 
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Electronic  Structure  of  the  Lowest 1  Aj  and  3Bi  States 
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Contribution  from  the  Department  of  Chemistry.  University  of  Utah. 
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Abatract:  Mu; ^configuration  telNconsisiem  field  calculations  have  been  performed  at  several  geometries  on  the  lowest  'Ai 
and  3B|  states  of  cyclopropenylidene.  The  effects  of  hybridization  and  electron  correlation  are  analyzed  at  the  various  geome¬ 
tries.  These  model  prototype  calculations  reveal  strong  correlation  effects  in  the  aromatic  singlet  state  as  well  as  a  singlet-trip¬ 
let  potential  surface  crossing  which  may  play  a  role  in  the  chemical  reactivity  of  analogous  chemical  reactions  involving  inser¬ 
tion  of  singlet  or  triplet  reagents  into  triple  bonds.  It  is  demonstrated  that  mulUconfigurational  wave  functions  are  absolutely 
necessary  to  obtain  a  qualitative  description  of  the  two  spin  states  of  these  unsaturated  carbencs. 


1.  Introduction 

Considerable  theoretical  ami  experimental  interest  has 
been  devoted  to  the  chemistry  of  saturated  and  unsaturated 
carbcnes.  Theoretical  research  in  our  laboratory'  and  theo¬ 
retical  and  experimental  work  in  other  laboratories1'3'11  has 
focused  on  the  very  important  question  of  singlet-triplet  energy 
differences  of  isolated  carbcnes.  Indeed,  postulation  of  one-step 

*  Camdl*  and  H<*r>  DrayfM  Tsaabar-tabalar,  OsvM  f.  Gardner  XiaaarcS 
F«Bo»,  John  Simon  fi»w. 


stereospecific  addition  of  singlet  carbenes  and  two-step  non- 
stereospecific  addition  of  triplet  carbenes  to  olefins  has  played 
a  dominant  role  in  the  experimental  determination  of  the  spin 
state  of  carbenoid  intermediates.1  In  this  report,  we  attempt 
to  gain  a  fundamental  understanding  of  a  particular  carbene 
reaction  that  appears  to  involve  a  nonstereospecific  singlet 
addition  of  a  cyclic  carbene.4  Diphenylcyclopropenylidene 
appears  to  insert  into  both  dimethyl  maleate  and  dimethyl 
fumarate  to  produce  the  same  spiropentene,  thereby  suggesting 
a  triplet-like  mechanism.  However,  the  carbene  itself  is  gen- 
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Polarization  Green’s  Function  with  Multiconfiguration 
Self-Consistent-Field  Reference  States 

AJIT  BANERJEE,  JOHN  W.  KENNEY.  Ill,  AND  JACK  SIMONS* 

Department  of  Chemistry,. University  of  Utah,  Salt  Lake  City,  Utah  84112,  U.S.A. 


Abstracts 

The  polarization  Green's-function  formalism  in  the  superoperator  notation  of  Gosanski  and 
Lukman  is  re-derived  using  a  multiconfiguration  self-consistent-field  (mc-scf)  reference  state  to 
establish  the  superoperator  metric.  The  potential  advantages  of  employing  this  more  general  reference 
state  in  Green's-function  theories  and  certain  inherent  weaknesses  associated  with  the  traditional 
Hartree-Fock  or  Rayleigh-Schrfidinger  perturbation  theory  reference  state  choices  are  briefly 
discussed.  The  Hermiticity  of  the  superoperators  is  analyzed  within  the  framework  of  the  mc-scf 
reference  state.  Using  a  concept  of  order  appropriate  for  this  reference  state  choice,  explicit  formulas 
and  computational  procedures  for  the  implementation  of  this  Green's-function  theory  are  presented 
and  specialized  to  include  terms  consistent  through  second  order. 

Le  formalisme  des  fonctions  de  Green  de  polarisation  dans  la  notation  utilisant  les  superopdrateurs 
de  Goscinski  et  Lukman  a  iti  rddtabli  avec  un  dtat  de  rdfdrence  mc-scf  pour  obtenir  le  mdtrique  des 
superopdrateurs.  Les  avantages  potentials  de  cet  itat  de  rdfdrence  plus  gdndnl  ainsi  que  certains 
points  faibles  inherent*  associds  aux  dtats  de  rdfdrence  utilisds  traditionnellement  dans  les  thdories  de 
perturbation  Hartree-Fock  ou  Rayleigh-Schrodinger  sont  discutds  bridvement.  L'hermhicitd  des 
superopdrateun  est  analysde  dans  le  cadre  des  dtats  de  rdfdrence  MC-SCF.  Utilisant  une  notion  d'ordrc 
propre  i  cet  dtat  de  rdfdrence  on  prdsente  des  formules  explicites  et  des  procddds  de  calcul  pour  la 
rdalisation  de  cette  thdorie  des  fonctions  de  Green,  qui  sont  spdcislisds  pour  indure  des  termet 
consistants  jusqu'au  second  ordre. 

Der  Formalismus  fur  Green’sche  Funktionen  von  Polarisierungstyp  in  der  Superoperator - 
bezeichnung  von  Goscinski  und  Lukman  wird  mittels  eines  MC-SCF-Referenzzustands  hergeleitet. 
Die  moglichen  Vorteile  mit  diesem  allgemeineren  Referenzzustand  und  gewisse  mit  den  in  den 
traditionellen  Hartree-Fock-  Oder  Rayleigh-Schrodinger-St&rungstheorien  verwendeten  Referenz- 
zustanden  verbundene  eigene  Schwichen  werden  kurz  diskutiert.  Die  Hermitizitit  der  Supero- 
peratoren  wird  im  Rahmen  des  MC-SCF-Referenzzustands  analysiert.  Mittels  eines  fur  diesen 
Referenzzustand  geeigneten  Ordnungsbegriffs  werden  explizite  Formeln  und  Rechnungsverfahren 
fur  die  Implementierung  dieser  Theorie  von  Green'schen  Funktionen  vorgelegt  und  spezialisiert  um 
Glieder  bis  zur  zweiten  Ordnung  einzuschliessen. 


I.  Introduction 

In  the  last  decade,  Green’s-function  approaches  to  the  direct  calculation  of 
electronic  excitation  energies  and  other  spectral  properties  have  been  employed 
with  considerable  success  to  many  atomic  and  molecular  systems  [1-18],  Electron 
attachment  and  detachment  processes  (i.e.,  ionization  processes)  also  have  been 
studied  extensively  with  Green’s-function  methods  [9-11,  19-27],  This  paper 
focuses  on  the  development  of  a  new  Green’s-function  approach  in  which  the 
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Analysis  of  the  “charge  resonance”  transition  in  anthracene 
dimer  anion 

Jeanne  McHale  and  Jack  Simons** 

Chemistry  Department.  University  of  Utah,  Salt  Lake  City.  Utah  84112 
(Received  20  August  1979;  accepted  24  September  1979) 

The  experimentally  observed  “charge  resonance”  transition  of  anthracene  dimer  anion  is  explained  in 
terms  of  the  resonance  interaction  of  the  lowest  unoccupied  orbitals  of  neighboring  anthracene  molecules. 

The  distortion  energies  of  the  anthracene  molecule  and  anion  are  shown  to  make  an  important 
contribution  to  the  excitation  energy.  The  numerical  estimate  of  the  resonance  integral,  based  od  the 
experimental  absorption  energy  of  the  dimer  anion,  is  shown  to  be  related  to  that  which  is  appropriate 
for  a  calculation  of  electron  mobility  in  crystalline  anthracene. 


I.  INTRODUCTION 

Shida  and  Iwata1  have  studied  the  electronic  absorp¬ 
tion  spectra  of  y-irradiated  solutions  of  dianthracene  in 
a  low -temperature  (77  ®K)  matrix  of  methyltetrahydro- 
furan  (MTHF).  Dianthracene  is  a  composite  molecule 
made  up  of  two  anthracene  molecules  in  a  sandwich  con¬ 
figuration,  with  a  bonds  connecting  the  equivalent  9  and 
10  positions  of  each  molecule.  Shida  and  Iwata  bleached 
dilute  y-irradiated  solutions  of  dianthracene  and  ob¬ 
served  the  resulting  optical  absorption  spectrum  to  be 
very  nearly  a  superposition  of  the  spectra  of  the  isolated 
anthracene  molecule  and  anthracene  anion.  In  addition, 
a  broad  absorption  occurring  between  about  5000  and 
10000  cm*1  (-0.6-1. 2  eV)  with  a  peak  at  6000  cm*1 
(0. 72  eV)  was  observed.  Since  the  absorption  in  this 
region  was  unlike  the  doublet  that  would  have  been  ob¬ 
served  for  dianthracene  anion,  and  because  the  rest  of 
the  spectrum  showed  peaks  which  could  be  assigned  to 
anthracene  molecules  and  anions,  Shida  and  Iwata  ex¬ 
plained  the  absorption  commencing  at  5000  cm*1  in  terms 
of  a  “charge  resonance"  transition  of  a  weakly  interact¬ 
ing  anthracene  dimer  anion  involving  an  anthracene 
molecule  and  anion  which  are  no  longer  held  together  by 
a  bonds.  The  dimer  anion  results  from  the  reaction  of 
mobile  electrons  (formed  by  bleaching  at  520  nm)  with 
the  o-bonded  dianthracene  anion.  According  to  Shida 
and  Iwata,  the  charge  resonance  (CR)  absorption  in¬ 
volves  as  upper  and  lower  states  the  symmetric  and  anti¬ 
symmetric  dimer  anion  wave  functions,  which  are  con¬ 
structed  from  the  ground-state  molecule  (0O)  and  anion 
(<£„)  wave  functions  at  sites  1  and  2: 

(i) 

Shida  and  Iwata  calculated  the  Intensities  and  positions 
of  peaks  in  the  spectrum  of  the  dimer  anion  using  the 
self -consistent -field  molecular  orbitals  (SCF-MO’s)  of 
anthracene  molecule  as  basis  orbitals,  assuming  the 
short  and  long  symmetry  axes  of  the  two  molecules  to 
be  parallel,  and  estimating  the  intermolecular  distance 
to  be  3.4  A  or  about  twice  the  van  der  Waals’  radius 
of  a  typical  s  system.  They  predict  that  the  charge 
resonance  absorption  will  occur  at  about  4000  cm*1, 
compared  to  the  experimentally  observed  energy  of 
5000  to  10  000  cm*1. 


^Camille  and  Henry*  Dreyfus  Fellow,  John  Simon  Guggenheim 
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In  Sec.  n,  we  show  how  the  CR  energy  can  be  more 
reasonably  accounted  for  in  terms  of  the  resonance 
energy  of  the  lowest  unoccupied  molecular  orbitals 
(LUMO’s)  on  the  neighboring  anthracenes  as  well  as  in 
terms  of  the  geometrical  distortion  energies  of  the 
planar  anion  and  nonplanar  molecule.  This  distortion 
energy,  which  shall  be  shown  to  make  an  important  con¬ 
tribution  to  the  CR  energy,  was  not  accounted  for  in  the 
calculation  of  Shida  and  Iwata,  who  assumed  planar 
geometries  for  both  molecules  in  the  dimer. 

II.  ANALYSIS  OF  THE  CHARGE  RESONANCE 
ENERGY  Af  CR 

In  the  following,  we  consider  the  CR  transition  to  in¬ 
volve  dimer  states  which  can  be  described  in  terms  of 
products  of  a  ground-state  anthracene  molecule  wave 
function  tut (r,  Q)  and  a  ground-state  anthracene  anion 
wave  function  Q).  The  subscripts  refer  to  sites  1 
and  2,  which  are  assumed  to  be  3. 4  A  apart.  To  consid¬ 
er  the  energy  of  the  vertical  CR  transition,  we  must 
keep  in  mind  that  the  equilibrium  geometry,  denoted  by 
the  internal  coordinates  Q,  is  not  necessarily  the  same 
for  the  anion  as  for  the  molecule.  Let  the  geometry 
which  minimizes  the  energy  of  the  molecule  be  repre¬ 
sented  by  Qp  (planar),  and  that  which  minimizes  the 
energy  of  the  anion  be  given  by  Qn  (nonplanar).  We  con¬ 
sider  the  zeroth-order  ground  state  to  consist  of  a 
planar  molecule  at  site  1  and  a  nonplanar  anion  at  site 
2: 

*l9}=*Ur,Q,)1>Hr,Qn).  (2) 

The  zeroth-order  excited  state  is  composed  of  a  planar 
anion  at  site  1  and  a  nonplanar  molecule  at  site  2: 

*r-«(r.Q#)^5(r,<?,).  (3) 

If  these  two  states  are  allowed  to  interact,  the  perturbed 
ground-  and  excited-state  energies  are  obtained  by 
diagonalizing  the  configuration  interaction  (Cl)  matrix 
H, 

To  simplify  the  calculation,  zero  overlap  of  wave  func¬ 
tions  on  different  sites  is  assumed.  In  this  approxima¬ 
tion,  the  Hu  matrix  element  is  Just  the  sum  of  the 
ground-state  equilibrium  energies  of  an  undistorted 
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Should  one  use  complex  basis  functions  in  coordinate 
rotation  calculations  on  molecules? 

Jack  Simons*' 
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In  applying  the  coordinate  rotation  method1  (CRM)  to 
electron-molecule1  scattering  resonances,  one  attempts 
to  find  the  complex  eigenvalues  (Z„)  corresponding  to 
square  integrable  eigenfunctions  of  the  rotated 
Hamiltonian  H„.  H„  is  obtained  from  the  conventional 
Hamiltonian  by  scaling  the  electronic  coordinates  r,  by 
the  complex  quantity  q(r,  -r,q).  For  a  diatomic  mole¬ 
cule  whose  nuclei  have  charges  ZA,  ZB  and  positions  RA, 
Rs.  contains  kinetic  energy  ( — tj ’/2  T,v  vrj),  electron- 
electron  interactions  (rjXuj  r',)),  and  electron-nuclei  in¬ 
teractions  l-vlt  ZAlr,  -q"'R a\'x  -nit  Zb\t,  -T)*1Ra!*1). 
The  desired  resonance  eigenvalues  of  Hr  are  usually 
found  by  expanding  i'„(r)  in  some  basis  of  square  inte¬ 
gral  functions  (t!„  aI,C»(q)<%(r))  which  are  themselves 
real  and  q -independent.  If  the  basis  set  {<5,(r)}  is  com¬ 
plete.  this  is  a  perfectly  reasonable  expansion  to  make. 
However,  in  practical  applications  of  the  CRM  to  elec  - 
tron-molecule  scattering,  the  set  {6*}  usually  consists 
of  a  finite  number  of  (A'  +  l)  electron  Slater  determinants 
involving  molecular  orbitals  expressed  in  terms  of 
atomic  orbitals  centered  at  R  .  and  Ra.  This  set  is  cer¬ 
tainly  not  complete,  so  it  is  quite  natural  to  ask  whether 
it  is  even  reasonable  to  employ  this  particular  choice  of 
basis . 

It  is  well  known’  that  the  lowest  energy  eigenfunctions 
of  the  unsealed  Hamiltonian  (/f(q  =  l)]  possess  high  am¬ 
plitude  near  each  of  the  two  nuclei.  This  arises  because 
the  Schrodinger  equation  is  dominated,  for  r,  =RA,  by 
_  r\  -  Z A I r,  -  RA I  *l  which,  just  as  in  the  hydrogen 
atom,  gives  large  (Is)  amplitude  in  r,  space  near  RA. 

In  choosing  atomic  orbital  basis  sets  for  use  in  standard 
variational  calculations  on  ff(q  =  1),  it  is  important  to 
include  functions  which  properly  represent  the  locally 
high  amplitudes  near  R  because  of  the  strong  influence 
of  this  part  of  the  electron  density  on  the  total  electronic 
energy.  It  is  most  common  to  employ  either  Is  Slater  or 
"tight"  contracted  Gaussian  basis  functions  centered  at 
R„  and  R(. 

In  considering  the  best  choice  of  basis  for  treating  the 
rotated  H„,  it  therefore  seems  important  to  guarantee 
that  the  resultant  approximations  to  ;'»(r)  display  high 


orbital  amplitude  at  points  dictated  by  B„.  Because  the 
electron-nuclei  interactions  in  Hn  have  singularities  at 
r,  =r)*lRx,a.  it  is  probably  quite  important  to  include  in 
any  reasonable  atomic  orbital  basis  (Slater  or  Gaussian) 
functions  centered  at  the  complex  positions  ij~1Ra  and 
n* ’Ra-  This  is  easily  achieved  in  electron-atom  prob¬ 
lems  by  simply  choosing  the  coordinate  origin  to  be  at 
the  atomic  nucleus;  then  q^R.,  =  R„  =  0.  For  molecules, 
this  choice  of  origin  does  not  remove  the  problem. 

Because  the  most  common  procedures1*1  for  locating 
resonance  energies  (Za)  within  the  CRM  involve  follow¬ 
ing  eigenenergies  of  If,  as  q  is  varied,  all  q  -dependent 
atomic -basis  integrals  (one-  and  two-electron)  arising 
in  the  computation  (e.g.,  by  configuration  interaction) 
of  Ztt  need  to  be  recomputed.  We  therefore  suggest  that 
one  employ  only  a  few  (e.g. ,  one  per  nucleus)  “tight” 
basis  functions  centered  at  q*1  RA,  q^Rg,  while  using 
large  numbers  of  atomic  orbitals  centered  at  RA,  R(. 

The  latter  basis  functions  should  be  capable  of  describ¬ 
ing  the  resonance  wave  function  ;-„(r)  in  the  “valence” 
region  and  for  values  of  r  which  are  on  the  "exterior” 
of  the  target  molecule.  Moreover,  only  those  integrals 
which  involve  the  functions  centered  at  q*JRA  „  B  need 
to  be  recomputed  as  q  is  varied. 

If  the  basis  {o,}  is  not  able  to  describe  high  orbital 
amplitudes  at  q'*Raia  then  spurious  imaginary  contribu¬ 
tions1  can  arise  in  the  desired  resonance  energies  and 
even  in  the  bound  states  which  should  be  q-independent. 
We  believe  that  the  "compromise  basis”  described 
above  has  potential  for  eliminating  the  undesirable  spu¬ 
rious  features.  Resclgno  et  al.*  and  Junker  et  al.* 
have,  in  fact,  already  employed  an  approach  which  is 
similar  to  that  being  proposed  here.  For  example, 
Rescigno  et  al.  employ  a  Slater  determinant  basis  in 
which  the  N  "bound”  orbitals  of  the  target  atom  are  not 
even  rotated  (in  our  language,  this  would  be  analogous  to 
expanding  the  A’  bound  orbitals  in  terms  of  the  basis  or¬ 
bitals  centered  at  rflRAtM,  which  essentially  "undoes” 
the  effort  of  rotating  these  A*  electrons  in  H,).  Only  the 
"active”  .V  +  1st  orbital  is  rotated.  There  is  only  a 
slight  difference  between  what  we  are  proposing  and  what 


Application  of  the  Coordinate  Rotation  Method 
to  Metastable  Atom-Diatom  Scattering 
Resonances 
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Abstract 

The  procedure  for  applying  the  coordinate  rotation  device  to  the  coupled-channel  scattering 
equations  is  presented.  Results  of  applications  of  the  method  to  resonances  occurring  in  low-energy 
atom-diatom  collision  processes  are  presented.  In  particular,  rotational  Feshbach  resonances  in 
model  (homogeneous)  potentials  which  have  also  been  studied  by  Crabenstetter  and  Le  Roy  are 
examined  within  this  novel  approach. 


1.  Introduction 

The  so-called  coordinate  rotation  method  (CRM)  (1]  has  been  used  success¬ 
fully  to  study  scattering  resonances  in  electron-atom  and  electron-molecule 
collisions  [2,  3].  More  recently,  it  has  also  been  used  to  study  rotational  pre¬ 
dissociation  in  HJ  [4].  In  the  present  work,  we  report  the  first  application  of  this 
technique  to  resonances  which  occur  in  closed-shell  atom-diatom  collisions. 
Unlike  the  HJ  rotational  predissociation  problem,  the  three-atom  situation 
addressed  here  involves  motion  on  a  three-dimensional  potential  energy  surface 
(HJ  involves  only  a  one-dimensional  curve).  As  a  result,  the  diatom’s  internal 
(vibrational-rotational)  energy  and  (rotational)  angular  momentum  must  be 
explicitly  coupled  to  the  corresponding  quantities  of  the  incident  atom.  The 
proper  coupling  of  energy  and  angular  momentum  (asymptotic)  basis  functions 
yields  the  usual  set  of  coupled-channel  equations.  These  coupled  differential 
equations,  upon  application  of  the  coordinate  rotation  device,  give  rise  to  a  set 
of  coupled-complex  differential  equations  for  the  components  of  the  scattering 
wavefunction  which  corres^iwd  to  various  channels. 

In  Sec.  2  we  obtain  the  set  of  coupled-channel  differential  equations  for  a 
closed-shell  atom-diatom  collision  and  we  demonstrate  how  these  equations 
are  transformed  upon  coordinate  rotation.  In  Sec.  3,  we  indicate  how  the  re¬ 
sulting  coupled-complex  differential  equations  can  be  solved  either  by  a  finite 
difference  method  or  by  an  expansion  of  the  (radial)  scattering  wave  amplitudes 
in  terms  of  square-integrable  basis  functions.  A  method  for  identifying  the  de¬ 
sired  metastable  resonance  energy  (which  is  complex)  and  its  accompanying 
wavefunction  is  also  outlined  in  Sec.  3.  Section  4  contains  our  concluding  re¬ 
marks. 
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Abstract 

The  procedure  for  applying  the  coordinate  rotation  device  to  the  coupled-channel  scattering 
equations  is  presented.  Results  of  applications  of  the  method  to  resonances  occurring  in  low-energy 
atom-diatom  collision  processes  are  presented.  In  particular,  rotational  Feshbach  resonances  in 
model  (homogeneous)  potentials  which  have  also  been  studied  by  Crabenstetter  and  Le  Roy  are 
examined  within  this  novel  approach. 


1.  Introduction 

The  so-called  coordinate  rotation  method  (CRM)  [1]  has  been  used  success¬ 
fully  to  study  scattering  resonances  in  electron-atom  and  electron-molecule 
collisions  [2,  3].  More  recently,  it  has  also  been  used  to  study  rotational  pre¬ 
dissociation  in  HJ  [4].  In  the  present  work,  we  report  the  first  application  of  this 
technique  to  resonances  which  occur  in  closed-shell  atom-diatom  collisions. 
Unlike  the  HJ  rotational  predissociation  problem,  the  three-atom  situation 
addressed  here  involves  motion  on  a  three-dimensional  potential  energy  surface 
(HJ  involves  only  a  one-dimensional  curve).  As  a  result,  the  diatom’s  internal 
(vibrational-rotational)  energy  and  (rotational)  angular  momentum  must  be 
explicitly  coupled  to  the  corresponding  quantities  of  the  incident  atom.  The 
proper  coupling  of  energy  and  angular  momentum  (asymptotic)  basis  functions 
yields  the  usual  set  of  coupled-channel  equations.  These  coupled  differential 
equations,  upon  application  of  the  coordinate  rotation  device,  give  rise  to  a  set 
of  coupled-complex  differential  equations  for  the  components  of  the  scattering 
wavefunction  which  correspond  to  various  channels. 

In  Sec.  2  we  obtain  the  set  of  coupled-channel  differential  equations  for  a 
closed-shell  atom-diatom  collision  and  we  demonstrate  how  these  equations 
are  transformed  upon  coordinate  rotation.  In  Sec.  3,  we  indicate  how  the  re¬ 
sulting  coupled-complex  differential  equations  can  be  solved  either  by  a  finite 
difference  method  or  by  an  expansion  of  the  (radial)  scattering  wave  amplitudes 
in  terms  of  square-integrable  basis  functions.  A  method  for  identifying  the  de¬ 
sired  metastable  resonance  energy  (which  is  complex)  and  its  accompanying 
wavefunction  is  also  outlined  in  Sec.  3.  Section  4  contains  our  concluding  re¬ 
marks. 

*  .University  of  Utah  Graduate  Research  Fellow. 

*  Camille  and  Henry  Dreyfus  Fellow  and  John  Simon  Guggenheim  Fellow. 


International  Journal  of  Quantum  Chemistry:  Quantum  Chemistry  Symposium  14. 467-475  (I9S0) 

?  19*0 by  John  Wiley  k  Sons.  Ine.  OI61-3642/SO/OOI4.0467SOI.OO 


The  Electronic  Structure  of  Singlet 
Cyclopentadienylidene 


RON  SHEPARD*  AND  JACK  SIMONS* 

Chemistry  Department,  University  of  Utah,  Salt  Lake  City,  Utah  84112,  US.  A. 

Abstract 

We  have  carried  out  ab  initio  mulliconfigurational  self-consistent  field  calculations  on  the  lowest 
singlet  state  of  the  carbene  C5H.4.  We  find  that  there  are,  along  a  restricted  geometry  deformation 
coordinate,  two  local  energy  minima  corresponding  to  localized  and  delocalized  x  bonds  having 
4r  and  6a-  electrons,  respectively.  The  role  of  4n  +  2  resonance  in  determining  the  nature  of  the 
electronic  structure  and  reactivity  of  this  carbene  is  analyzed  and  compared  with  that  of  other 
carbenes. 


1.  Introduction 

In  previous  studies  of  the  electronic  structure  of  linear  and  cyclic  carbenes 
[1-3],  the  importance  of  orbital  hybridization  of  the  electron-deficient  carbene 
carbon,  electron  delocalization,  and  electron  correlation  have  been  investigated 
in  order  to  gain  a  basic  understanding  of  the  nature  of  these  important  reaction 
intermediates.  The  most  important  correlation  configuration  of  the  XA  \  states 
(within  Cy  geometries)  of  methylene  (CH2),  cyclopropenylidene  (C3H2),  and 
linear  unsaturated  carbenes  can  be  described  as  a  double  excitation  from  an 
occupied  orbital  into  a  vacant  b\  orbital  of  the  dominant  configuration.  A 
comparison  [1  ]  of  the  importance  of  this  a?  -»■  b]  configuration  in  the  cyclic  C3H2 
and  the  other  systems  [2, 3]  shows  that  this  correlation  is  less  important  in  C3H2. 
In  the  C3H2  system,  the  b]  configuration  has  4 it  electrons  and  thus  acquires 
antiaromatic  character  while  the  dominant  a]  configuration  has  2zr  electrons 
and  hence  possesses  aromatic  stability  because  it  satisfies  the  “4 n  +  2  rule.” 
This  indicates  that  for  this  cyclic  carbene  the  aromatic  nature  of  the  correlation 
configuration  must  be  considered  along  with  the  effects  of  delocation  and  hy¬ 
bridization. 

In  this  brief  report,  the  importance  of  aromatic  factors  is  investigated  for 
another  cyclic  carbene — cyclopentadienylidene.  In  contrast  to  the  C3H2  carbene, 
both  the  delocalization  effects  and  the  aromatic  nature  of  the  aj  config¬ 
uration  work  in  accord  to  increase  the  importance  of  this  configuration.  The 
result  is  that  the  b 7  configuration  becomes  not  only  more  important  than  in  CH2 
and  C3H;  but  it  actually  becomes  the  dominant  configuration  at  some  molecular 
geometries  examined  in  this  study. 

*  Present  address:  Bauellc  Memorial  Institute.  Columbus.  Ohio  43201. 
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Abstract 

We  have  successfully  used  graphical  unitary  group  methods  to  implement  efficiently  both  mul- 
ticonfiguralional  self-consistent  field  and  configuration  interaction  ab  initio  quantum  chemical 
computer  programs.  We  indicate  how  these  group  theoretic  methods  improve  efficiency  and  we 
demonstrate  their  utility  on  the  ground  and  excited  states  of  a  model  problem  for  which  a  full- 
configuration  interaction  calculation  can  be  performed. 


I.  Introduction 

Using  the  exponential-Mambda  (ElL)  method  [1]  for  wavefunction  optimi¬ 
zation  and  the  graphical  unitary  group  method  [2]  for  the  calculation  of  one- 
and  two-particle  density  matrix  elements,  we  have  performed  multiconfiguration 
self-cons.stent  field  (MCSCF)  calculations  on  the  ground  and  excited  states  of 
the  Be  ■<  mat  several  levels  of  sophistication  as  a  model  calculation.  We  show 
how  the  global  convergence  problems  of  the  ElL  method  may  be  overcome  in 
these  cases.  We  further  show  that  although  completely  flexible  configuration 
selection  is  not  easily  achieved  within  the  unitary  group  formalism  which  we 
have  implemented,  limited  distinct  row  tables  based  only  on  occupation  re¬ 
strictions  are  easily  constructed.  The  resulting  configuration  lists  are  neither 
so  large  that  chemical  insights  are  lost  nor  are  they  so  restrictive  that  chemical 
accuracy  is  sacrificed.  We  show  how  the  graphical  unitary  group  method  may 
be  used  to  construct  efficiently  the  matrix  elements  required  in  the  ElL 
MCSCF  method  without  explicitly  performing  the  two-particle  density  matrix 
construction.  Because  the  ElL  method  possesses  very  good  local  convergence 
properties,  not  only  accurate  energy  expectation  values  but  also  very  accurate 
MCSCF  wavefunctions  are  reported. 

The  MCSCF  wavefunction  is  a  multiconfiguration  wavefunction: 

|*mc>  »  £C*|tf>  (1) 

A 

where  both  the  configuration  mixing  coefficients  C *  and  the  orbitals  are  opti¬ 
mized  to  give  the  lowest  energy  expectation  value.  A  property  of  the  MCSCF 
wavefunction  is  that  it  satisfies  the  generalized  Brillouin  theorem  [3]  (GBT): 
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Abstract* 

The  coupled-duster  approach  to  obtaining  the  bond-state  wave  functions  of  many-electron 
systems  is  extended,  with  a  set  of  physically  reasonable  approximations,  to  admit  a  multiconfiguration 
reference  state.  This  extension  permits  electronic  structure  calculations  to  be  performed  on  correlated 
closed-  or  open-shell  systems  with  potentially  uniform  precision  for  all  molecular  geometries.  Explicit 
coupled  cluster  working  equations  are  derived  using  a  multiconfiguration  reference  state  for  the  case  in 
which  the  so-called  cluster  operator  is  approximated  by  its -one-  and  two-particle  components.  The 
evaluation  of  the  requisite  matrix  elements  is  facilitated  by  use  of  the  unitary  group  generators  which 
have  recently  received  wide  attention  and  use  in  the  quantum  chemistry  community. 

Le  procede  des  amas  couples  pour  obtenir  des  fonctions  d'onde  pour  les  itats  lies  d'un  systime  i  N 
Electrons  a  it6  generalise  dans  le  cadre  d’un  nombre  d'approximations  raisonnables  du  point  de  vue 
physique,  pour  utiliser  un  etat  de  rifitence  multiconfigurationnel.  Cette  extension  permet  des  calculi 
de  la  structure  electronique  pour  des  systimes  i  couches  fermies  ou  ouvertes  correldes  avec  une 
precision  qui  est  uniforme  pour  toutes  les  geometries  moldculaires.  Des  equations  explicates  sont 
obtenues  avec  un  etat  de  reference  multiconfigurationnel  dans  le  cas  ou  le  soi-disant  opdrateur  d’amas 
est  remplace  par  ses  composantes  1  une  et  deux  particules.  Le  calcul  des  elements  matriciels 
ndcessaires  est  facilite  par  I'utilisation  des  gendrateurs  du  groupe  unitairc. 

Das  Verfahren  mit  gekoppelten  “Clusters"  fur  die  Berechnung  von  Wellenfunktionen  fur  gebun- 
dene  Zustande  eines  Vielelektronensystems  ist  im  Rahmen  gewisser  physikalisch  angemessenen 
Naherungen  erweitert  worden,  um  einen  multikonfigurationellen  Referenzzustand  zu  erlauben.  Mit 
dieser  Erweiterung  konnen  Berechnungen  der  Elektronenstruktur  von  Systemen  mit  korrelierten 
abgeschlossenen  oder  offenen  Schalen  und  mit  einer  potentiell  einheitlichen  Genauigkeit  fur  alle 
Molekiilgeometrien  ausgefuhrt  werden.  Explizite  Gleichungen  warden  mit  einem  multi¬ 
konfigurationellen  Referenzzustand  fur  den  Fall  abgeleitet.  wo  der  sogenannte  Gusteroperator  von 
seinem  Ein-  und  Zweiteilchenkomponenten  angenahert  wird.  Die  Auswertung  der  erforderlichen 
Matrixelcmente  wird  durch  die  Anwendung  der  Generatoren  der  unitiren  Gruppe  erleichtert. 


1.  Introduction 

A  cluster  expansion  formalism  was  developed  by  Coester  [1]  for  treating 
S-  fermion  nuclear  systems.  Cizek  [2]  and  others  [3]  expressed  this  formalism 
such  that  it  was  more  suitable  for  electronic  structure  calculations  on  closed-shell 
systems  and  actually  applied  the  resulting  equations  in  an  ab  initio  framework  to 
the  ground  states  of  a  few  closed-shell  molecules  [4).  More  recently,  attempts 
have  been  made  to  extend  the  formalism  to  open-shell  systems  [5, 6]  and  Freeman 
et  al.  [7]  have  applied  the  method  to  the  uniform  electron  gas. 
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